Andreev Transport through Side-Coupled Double Quantum Dots 
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We study the transport through side-coupled double quantum dots, connected to normal and 
superconducting (SC) leads with a T-shape configuration. We find, using the numerical renormal- 
ization group, that the Coulomb interaction suppresses SC interference in the side dot, and enhances 
the conductance substantially in the Kondo regime. This behavior stands in total contrast to a wide 
Kondo valley seen in the normal transport. The SC proximity penetrating into the interfacial dot 
pushes the Kondo clouds, which screens the local moment in the side dot, towards the normal lead 
to make the singlet bond long. The conductance shows a peak of unitary limit as the cloud expands. 
Furthermore, two separate Fano structures appear in the gate-voltage dependence of the Andreev 
transport, where a single reduced plateau appears in the normal transport. 

PACS numbers: 73.63.Kv, 74.45.+C, 72.15.Qm 
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I. INTRODUCTION 

Observation of the Kondo effect in a quantum dot 
(QD)i has stimulated researches in the field of quantum 
transport, and recent experimental developments enable 
one to examine the Kondo physics in a variety of sys- 
tems, such as an Aharonov-Bohm (AB) ring with a QD 
and double quantum dots (DQD). In these systems mul- 
tiple paths for electron propagation also affect the tun- 
neling currents, and the interference causes Fano-type 
asymmetric line shapes. 

Superconductivity also brings rich and interesting 
features into the quantum transport. Competition 
between superconductivity and the Kondo effect has 
been reported to be observed in carbon nanotube 
QD and in semiconductor nanowires^&^&i Further- 
more, interplay between the Andreev scattering and 
the Kondo effect has been studied intensively for a 
QD coupled to a normal (N) lead and supercon- 
ductor (S), theoretically^^^^^^^ and 
experimentally^ So far, however, the Andreev-Kondo 
physics has been discussed mainly for a single dot . In this 
paper, we consider a DQD system with a T-shape geom- 
etry as shown in Fig. [T] and study how multiple paths 
affect the interplay at low temperatures, using the numer- 
ical renormalization group (NRG) method^ Golub and 
Avishai calculated first, to our knowledge, the Andreev 
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FIG. 1: (color online) Quantum dots coupled to (a) a normal 
lead (N) and superconductor (S), and (b) two normal leads. 
QD1 (QD2) is referred to as the interfacial (side) dot. Dashed 
line illustrates a dominant singlet pair in the Kondo cloud. 



transport through an AB ring with a QD^ in which a 
similar interference effect is expected. However, the un- 
derlying Kondo physics in such a combination with super- 
conductivity and interference is still not fully understood, 
and is needed to be clarified precisely, as measurements 
are being not impossible^ 

We find that the Coulomb interaction in the side dot 
(QD2 in Fig. [T]) suppresses destructive interference typ- 
ical of the T-shape geometry, and it enhances substan- 
tially the tunneling current between the normal and su- 
perconducting (SC) leads in the Kondo regime. This 
is quite different from the behavior seen in the nor- 
mal transport in the same configuration Fig. [1] (b), 
for which the conductance is suppressed, and shows 
a wide minimum called a Kondo valley as a result 
of strong interference by the Kondo resonance in the 
side dotM^2^^L2^1 The SC proximity penetrat- 
ing into the interfacial dot (QD1 in Fig. Q} causes this 
stark contrast between the Andreev and normal trans- 
ports. It also changes the Fano line shape in the gate- 
voltage dependence of the conductance. Furthermore, 
we show that the proximity deforms the Kondo cloud to 
make a singlet bond long, and it can be deduced from 
the Fermi-liquid properties of the ground state. 

In Sec. HH we introduce the model and describe the 
effective Hamiltonian in a large gap limit. In Sec. IIII) 
we show the numerical results and clarify the transport 
properties using the renormalized parameters. The Fano 
structures in the gate- voltage dependence of the conduc- 
tance are also discussed. A brief summary is given in the 
last section. 



II. MODEL 

We start with an Anderson impurity connected to SC 
and normal leads, 



H — Hdqd + Hs + Hn + Ht.s + Ht,n, 



(1) 



2 



where 

Hdqd = E { + y) - 1) + y - I) 2 } 

i=l,2 ^ ' ' 

+ 1 E (4^ + H - c -) > 

Hs = J2 £ k c SM c SM - E ( A 4,fcT C S,-H + H - c -) - 

fc,<7 

^=E^( C U^ + H.c), = (2) 

Hdqd describes the interfacial (i = 1) and side (i = 2) 
dots: £d,i the energy level, Ui the Coulomb interaction, 
n, = dl a d ia , and £ the inter-dot hopping matrix cl- 
ement. H s /n describes the SC/normal lead, and A is a 
s-wave BCS gap. Vs/n l& the tunneling matrix element 
between QD1 and the SC/normal lead. We assume that 
rgyjv(s) = ^s/n Sfe <5( £ ~ e k)/N is a constant indepen- 
dent of the energy e, where M is the total number of fc's 
in the leads. Throughout the work, we concentrate on a 
large gap limit A — > oo. Then the starting Hamiltonian 
H can be mapped exactly onto a single-channel model, 
which still captures the essential physics of the Andreev 
reflection and makes NRG approach efficient ; 19 ! 30 ' 31 

H cS = H s s + Hdqd + Ht,n + H^ , (3) 

Hf = - A dl (4Ai + H - c -) ' ^ 

A d i = r s . (5) 

Note that at A — -> oo the real and virtual excitations 
towards the continuum states outside the gap in the SC 
lead are prohibited. Nevertheless, the proximity from the 
SC lead to the dot remains finite, and it induces a local 
static pair potential A^i (= Ts) at QD1. Furthermore, 
the current can flow between the SC lead and the QD1 
via Ts- 



III. NUMERICAL RESULTS 
A. £ d2 -dependence 

We can calculate the conductance at zero temperature 
as a function of the level position Sd2 of QD2 for different 
values of U2, using the Kubo formula^ In this paper we 
focus on the Coulomb interaction in the side dot (QD2), 
assuming that U\ = in the following. The results of the 
conductance are shown in Fig. for e&\ = Ui = and 
Adi = l.Ot. The coupling to the normal lead is chosen to 
be (a) Tn — l.Ot and (b) Tjv = 0.2t. The conductance 
is enhanced for the Kondo regime — U2 < £<n < 0, where 
the wide Kondo valley appears in the case of the normal 
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FIG. 2: (color online) Conductance vs EdijU-i for (a) Fn = 
l.Ot and (b) Tjv = 0.2t, for several side-dot repulsions U2- 
The parameters for QD1 are chosen to be Edi =U% =0, and 
Adi = l.Ot which is the local SC gap defined by Adi = Ts- 



transport. This is a novel phenomenon caused by the in- 
terplay between superconductivity and the Kondo effect; 
the local gap A<ji due to the proximity into QD1 leads to 
the Andreev transport with destructive interference, but 
the introduction of U2 suppresses the SC interference via 
QD2, which in turn enhances the conductance. Note that 
the couplings are symmetric = for Fig. [2] (a), and 
in this particular case the conductance increases with U2 
for any Sd2, except for the values Ed2 — —U2 and 0.0. 
Outside of the Kondo regime, the side dot is empty or 
doubly occupied, and the interference becomes no longer 
important. When the coupling to the normal lead is small 
Tn < l.Ot as Fig. [2] (b), the conductance in the Kondo 
regime decreases after the peak reaches the unitary limit 
4e 2 //i. This behavior can be related to a crossover from 
short-range to long-range Kondo screening as illustrated 
in Fig. [H and is discussed later again. 

We examine the behavior at the middle point 6d2 = 
—U2/2 of Fig. [2] in detail. The low-lying energy states 
show the Fermi-liquid properties^ and the conductance 
can be deduced from the renormalized parameters for the 
quasi-particles (see Eq. ([7])). The conductance is plotted 
as a function of Tg/t in Fig. [3] for (a) Tn = l.Ot and 
(b) Tn = 0.2t. We see that the peak shifts towards 
smaller T$ as U2 increases and will coincide in the limit 
of U2 —> 00 with the dashed line, which corresponds to 
the conductance without the side dot. It means that the 
interference caused by the side dot is suppressed com- 
pletely for large U2, and in this limit the conductance 
reaches the unitary limit value for the symmetric cou- 
plings Ts — Tn. The difference in the line shape of 
Fig. [H (a) and that of Fig. H (b) at fixed T s reflects the 
position of the unitary-limit peak in Fig. [3] 



3 




FIG. 3: (color online) Conductance vs Ts/t at Ed2 = — U2/2 
and £di = Ui = for several values of U2. (a) Tat = l.Ot 
and (b) Tiv = 0.2t. The dashed line is the conductance of 
a single dot (t = 0), for which the horizontal axis should be 
interpreted as (a) 1.0r s /Tjv and (b) 0.2r s /riv. 




B. Fermi-liquid properties at Sd2 



-U2/2 



In order to clarify the properties of the the ground 
state precisely, we consider a special case Ed2 = —U 2 /2. 
Then the Hamiltonian H eS can be written in terms of the 
interacting Bogoliubov particles, which conserve the total 
charge, as shown in Appendix.— Consequently, the low- 
energy states can be described by a local Fermi liquid, 
the fixed-point Hamiltonian^ of which can be written in 
the form 



- Adi (4t4i +H.C.1 +H T 



H.c.) + t J2 ( d ld 2a + H.c 



(6) 



Here, Ad2 is a local SC gap that emerges in QD2 via the 
self-energy correction due to U2, while A^i = T$ as de- 
fined in Eq. ([5]) is caused by the direct proximity from the 
SC lead, t is the renormalized value of the inter-dot hop- 
ping matrix element. We calculate these parameters from 
the fixed-point of NRG.— Then, the conductance G and 
a staggered sum K of the pair correlation are deduced 
from the phase shift, 9, of the Bogoliubov particles, 



G 



4e 2 



tant 



sm z 29 



Ad2 Tjv 



K 



£(-1)^ 



2 = 1,2 



26» 

7T 



(7) 



t 2 - A d2 r s 



(4Ai 



^T>- ( 8 ) 



In Fig. HI we show the U2 dependence of the ground- 
state properties at £^2 = —U2/2 for e dl = U\ = and 
= l.Oi. The coupling is chosen to be Tjv = 0.2£, and 



FIG. 4: (color online) Ground state properties at e d2 = 
— U2/2: (a) Conductance, (b) t, A d2 , (c) «i, k 2 , and K = 
hi — K2- We choose sai = Ui = 0, Fg = l.Ot, and for filled 
(open) marks Vn = 1-Ot (0.2t). Inset of (b): Phase bound- 
ary, between singlet and doublet ground states, for an isolated 
DQD (rjv = 0) with a finite local SC gap A dl = T S - 



l.Ot. The conductance for IV < Ts shows a peak as 
a function of U 2 , while for Tn = Tg it increases simply 
towards the unitary limit. This corresponds to the differ- 
ence that we see in Fig.^a) and (b) at ed 2 = —U2/2. Fig- 
urclDJb) shows the renormalized parameters (o,«) t and 
(A, A) A 2 . The ratio t/t, which equals to the square root 
of the wavefunction renormalization factor Z (see Ap- 
pendix), decreases monotonically from 1.0 to 0.0 with in- 
creasing U2, while the local SC gap A d2 becomes large for 
intermediate values of U 2 . The behavior of these Fermi- 
liquid parameters implies that a crossover from weak to 
strong correlation regimes occurs around U 2 ~ 4.5£. The 
nature of the crossover can be related to a level crossing 
taking place in a molecule limit Tjy = 0, where QD1 is de- 
coupled from the normal lead. In this limit, the isolated 
DQD is described by a Hamiltonian Hbqd , , which 
includes the local SC gap A d i = Ts at QD1. The ground 
state of the molecule is a singlet or doublet depending on 
U2/t and Ts/t, as shown in the inset of Fig. [31(b) . The 
ground state is a spin-singlet, if either U 2 /t or Ts/t is 
small. In the opposite case, a spin-doublet becomes the 
ground state. The local moment in this doublet state 
emerges mainly at QD2, because the correlation at QD1 
is small in the present situation Sdi = U\ = 0. We see in 
the phase diagram in Fig. Sib) that the transition takes 
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place in this molecule limit at U2 — 4.5i for T$ = 1.04, 
and it agrees well with the position where the conduc- 
tance peak appears in Fig. Ufa). 

For finite Tjv, conduction electrons can tunnel from the 
normal lead to QD2 via QD1. However, the SC correla- 
tion Adi = Ts tends to make the local state at QD1 a sin- 
glet, which consists of a linear combination of the empty 
and doubly occupied states. Thus, for large Tg, the elec- 
trons at QD1 can not contribute to the screening of the 
moment at QD2. In this situation, the Kondo screening 
is achieved mainly by the conduction electrons tunnel- 
ing to QD2 virtually via QD1. This process is analogous 
to a supcrexchange mechanism, which can also be ex- 
pected from the Hamiltonian written in terms of the Bo- 
goliubov particles (|A.6|) in Appendix. From these obser- 
vations we see that the conductance peak at U2 — 4.5f in 
Fig. Ela) reflects the crossover from the short-range sin- 
glet to long-range one due to the superexchange screen- 
ing process (see also Fig. [1| for the Bogoliubov particles. 
Note that the peak structure of the conductance vanishes 
for T N = T s . 

The deformation of the Kondo cloud can also be de- 
duced from the results shown in Fig. H{c). This is be- 
cause the staggered pair correlation K is related directly 
to the scattering phase shift 9 of the Bogoliubov par- 
ticles, by the Fricdel sum rule Eq. (|A.10[) given in Ap- 
pendix. Therefore, the value of K reflects the changes 
occurring in the Kondo clouds. Particularly, a sudden 
change observed in K around U2 — 4.5£ shows that the 
phase of the wavefunction shifts by A9 ~ OAir during this 
change. This also explains the occurring of the crossover 
from the short-range to long-range screening. We can 
also calculate each correlation function m directly with 
NRG based on the definition. The local SC correlations 
K\ and k-2 have the same value in the nonintcracting case 
U2 = £d2 = 0, and thus in this particular limit there is no 
reduction in the amplitude of the proximity from QD1 to 
QD2. The Coulomb interaction U2 causes the reduction, 
as both m and K2 decrease for small U2 where the ground 
state is a singlet with a molecule character. For large U2, 
the SC correlation K2 almost vanishes in QD2, while k± 
shows an upturn and approaches the value expected for 
t = 0. Therefore, the SC proximity into QD1 is enhanced 
when the Kondo cloud expands to form a long-range sin- 
glet. Then the tunneling current is not interfered much 
by the local moment at QD2, and flows almost directly 
without using the path to the side dot. 



C. Fano line shape for Sdi 7^ 

So far, we have chosen the level of QD1 to be e^i = 0. 
The result obtained for different values of Edi is plotted 
as a function of Sd2 in Fig. [5] We see that two asymmet- 
ric Fano structures, each of which consists of a pair of 
peak and dip, emerge at Ed2 — and — U2 for Sdi 7^ 0, 
as the Fermi level crosses the energy corresponding to 
the upper and lower levels of the atomic limit. The con- 
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FIG. 5: (color online) Conductance vs Ed2 for several values 
of Edi, where U\ = and Ts = l.Ot. Other parameters: (a) 
Tjv = l.Oi, U 2 = 8.0*. (b) F N = 0.2*, U 2 = G.5t. 



ductance peaks become sharper for Tjv < Tj as |edi| 
increases. Maruyama et al studied the Fano structure in 
the side-coupled DQD system with two normal leads^ 
and showed that at low temperatures the conductance 
has only one pair of the peak and dip, which are sepa- 
rated widely by a Fano-Kondo plateau at — C/ 2 < £d2 < 0. 
This type of the plateau was known earlier to appear in 
a QD embedded in an AB ringi^ In contrast, our result 
shows that the Fano-Kondo plateau vanishes, when one of 
the leads is a superconductor. This is because the Kondo 
screening in this case is achieved by the long-range sin- 
glet bond due to the superexchange process, as a result 
of the competition between the SC proximity into QD1 
and Coulomb interaction at QD2. 



IV. SUMMARY 

We have studied Andreev transport through the side- 
coupled DQD with NRG approach. We have found that 
the Coulomb interaction in the side dot suppresses the 
destructive interference effect typical of the T-shape ge- 
ometry, and enhances the tunneling current between the 
normal and SC leads. This novel phenomenon is caused 
by the interplay between the SC correlation and the 
Kondo effect; the SC proximity into QD1 pushes the 
Kondo cloud towards the normal lead, and the conduc- 
tance shows a peak of the unitary limit as the nature of 
the singlet changes from a short-range to long-range one. 
We have also clarified that two asymmetric Fano struc- 
tures appear in the gate-voltage dependence of the An- 
dreev transport, instead of a reduced single Fano-Kondo 
plateau which appears in the Kondo regime of the normal 
transport. 
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APPENDIX: BOGOLIUBOV PARTICLES 

The Hamiltonian H eS defined in Eq. © can be trans- 
formed, at Ed2 = —U2/2, into the interacting Bogoliubov 
particles, which conserve the total charge. For describing 
this property briefly, we rewrite H using the logarith- 
mic discretization of NRGi^I 



Hf RG = K {N ~ 1)/2 (Hf + H DQD + H T> . 

N-l 

H.T,N + H-N = ^2 E *nA~ n/2 {fn+lafna 



hH.c). 



(Al) 



For n > 0, the operator f na describes the conduction 
electron in the normal lead, and t n is given by 



tn = D 



1 + 1/A 



1 - 1/A n+1 



y/1 - 1/A 2 »+Vl - 1/A 2 "+ 3 



(A.2) 



Here, D is the half- width of the conduction band. For 
the double-dot part, we use a notation /_j )CT = di a for 
i = 1, 2. Correspondingly, i_i = vA -1 / 2 and i_2 = t A -1 
with 



2T N DA A 



1 



VAN 



1-1/Ay 



log A. (A.3) 



At £d2 = £d2 + U2/2 = 0, the system has a uniaxial 
symmetry in the Nambu pseudo-spin spaced and the 
Hamiltonian can be simplified by the transformation 



(A.4) 
(A.5) 




Here, n 7i = Y^allaliai an< l tnc t°t a l number of the Bo- 
goliubov particles, Q 1 = 52iL-2^7«' * s conserved. The 
equation (|A.6[) implies that the low-energy excited states 
can be described by a local Fermi-liquid theory. This is 
true also for the original Hamiltonian H cS , and it does 
not depend on the discretization procedure of NRG.— 

To be specific, we assume that U\ = in the follow- 
ing. In this case, the Bogoliubov particles feel a nor- 
mal impurity potential E d .\ at QD1, and this poten- 
tial causes the superexchange mechanism that makes the 
singlet-bond long range as discussed in Scc. lIII Bl The re- 
tarded Green's function for the Bogoliubov particle 7-2,0- 
at QD2 takes the form 



G 7 {u) 



(A.7) 



lo — £ 7 (w) — 



lo - E dl + iT N 



where E 7 (w) is the self-energy caused by the interaction 
(U2/2) (^7,-2 — l) 2 - At zero temperature, the asymp- 
totic form of the Green's function for small ui ~ is 
given by 



Z 



V<£2 



U3-E. 



ill 



A d2 = ZS 7 (0), t = Vzt, z- l =\ 



9E 7 (w) 



(A. 



^=0 
(A.9) 



Note that A d 2 has a finite value even though ^2 = 0, 
because E d \ 7^ 0. The value of the parameters A^2 and 
t can be deduced from the fixed point of NRG^ Then, 
using the Friedel sum rule for Eq. (|A.6|) , the local charge 
at the double dot can be calculated from the phase shift 
6 of the Bogoliubov particles, 



"7,-2 



™7,-l 



9 = tan 



7T 

A rf2 r 



(n-6) 



(Z2 1 N 



A d2 E dl 



(A.10) 
(All) 



Here, E dx = y/^ + |A dl | 2 , £ dl = e dl + Ui/2, and 
A d i = Ts as defined in Eq. (|SJ). Then is trans- 

formed into a normal two-impurity Anderson model for 
the Bogoliubov particles 



WSFrg = A (JV - 



l)/2 



Edi (n 7 



1) 



E 

i=l,2 



N-l 

E 



E 



-n/2 



(7,1+1, 



H.c, 



-i) a 



(A.6) 



The charge of the Bogoliubov particles corresponds to the 
SC pair correlation for the original electrons f na . Specif- 
ically for £<0 = 0, it is transformed into the staggered 
sum K given in Eq. 0, by the inverse transformation of 
Eq. (|A.4[) . Similarly, the conductance G can be expressed 
in terms of the phase shift 8. Furthermore, the free quasi- 
particles corresponding to the Green's function given in 
Eq. (|A.8[) can be described by a Hamiltonian, which is 
rewritten in terms of the original electron operators in 
Eq. ([S]) by the inverse transformation. 
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